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Abstract 

Grobner bases can be used for computing the Hilbert basis of a 
numerical submonoid. By using these techniques, we provide an al- 
gorithm that calculates a basis of a subspace of a finite-dimensional 
vector space over a finite prime field given as a matrix kernel. 
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1 Introduction 

Grobner bases provide a uniform approach to tackling a wide range of prob- 
lems such as the solvability and solving algebraic systems of equations, ideal 
and radial membership decision, and effective computation in residue class 
rings modulo polynomial ideals [H [21 |6l [12] . 

Furthermore, Grobner basis techniques are not only a powerful tool for 
the algorithmic solution of some fundamental problems in commutative al- 
gebra they also provide means of solving a wide range of problems in 
integer programming and invariant theory once these problems have been 
expressed in terms of sets of multivariate polynomials |5l [TOl [13] . One such 
problem is the computation of the Hilbert basis for a submonoid of the nu- 
merical monoid Nq. This problem can be written in terms of polynomials 
and then be solved using Grobner basis techniques (TO]. Other elaborations 
of this method can be found in [3 [13] . 
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In this paper we will establish an algorithm using Grobner basis tech- 
niques that allows to calculate a basis for a subspace of a finite-dimensional 
vector space over a finite prime field given as a matrix kernel. This algorithm 
is based on the one for computing Hilbert bases proposed in [13] and is moti- 
vated by the fact that linear codes can be described as such subspaces [9l[T^. 

This paper is organized as follows. The first section provides an introduc- 
tion to Grobner bases, Hilbert bases and their construction for a submonoid 
of the numerical monoid Nq, and linear codes. The second section contains 
the main theorem and a variant of the algorithm for computing a basis for a 
subspace of described as a matrix kernel, where p is a prime. The paper 
concludes with an example illustrating the algorithm and its application to 
linear codes. 

2 Preliminaries 

Throughout this paper, Z denotes the ring of integers, Nq stands for the set of 
non-negative integers, IK denotes an arbitrary field, and ]K[x] = ]K[xi, . . . , x„] 
is the commutative polynomial ring in n indeterminates over K. 

2.1 Grobner bases 

The monomials in ]K[x] are denoted by x" and are identified 

with the lattice points u = (mi, . . . G Nq. The degree of a monomial 
x" is the sum |u| = ui + ■ ■ ■ + Un and the degree of a polynomial / is the 
maximal degree of all monomials involved in /. A term in K[x] is a scalar 
times a monomial. 

Denote by K[xf ^, . . . , xj^] the set of all polynomials given by monomials 
with exponents in Z", which is called the ring of Laurent polynomials. Nega- 
tive exponents can be overcome by introducing an additional indeterminate t. 
More precisely, we have 

K[xf\ Xf] ^ K[X1, ...,Xn,t]/ {XiX2 ...Xnt-1). (l) 

A monomial order on K[x] is a relation >- on the set of monomials x" in 
K[x] (or equivalently, on the exponent vectors in Ng) satisfying: (1) >- is a 
total ordering, (2) the zero vector is the unique minimal element, and (3) 
u >- V implies u + w >- v + w for all u, v, w G Nq . Familiar monomial orders 
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are the lexicographic order, the degree lexicographic order, and the degree 
reverse lexicographic order. 

Given a monomial order >-, each non-zero polynomial / G K[x] has a 
unique leading term, denoted by lt^(/) or simply lt(/), which is given by the 
largest involved term. The coefficient and the monomial of the leading term 
are called the leading coefficient and the leading monomial, respectively. 

If / is an ideal in ]K[x] and is a monomial order on ]K[x], its leading ideal 
is the monomial ideal generated by the leading monomials of its elements. 



A finite subset Q of an ideal / in K[x] is a Grobner basis for I with respect 
to >- if the leading ideal of / is generated by the set of leading monomials 
in Q; that is. 



If no monomial in this generating set is redundant, the Grobner basis will be 
called minimal. It is called reduced if for any two distinct elements g,h & Q, 
no term of h is divisible hy \i{g). A reduced Grobner basis is uniquely 
determined provided that the generators are monic. 

A Grobner basis for an ideal / in K[x] with respect to a monomial or- 
der >- on ]K[x] can be calculated by Buchberger's algorithm. It starts with 
an arbitrary generating set for I and provides in each step new elements of / 
yielding eventually a Grobner basis, which can further be transformed into a 
reduced one. For more about Grobner basics the reader may consult [IIEIE]- 

2.2 Monoids, Hilbert bases and their computation us- 
ing Grobner bases 

A monoid is a set M together with a binary operation such that the opera- 
tion is associative and M possesses an identitiy element. A submonoid of a 
monoid M is a subset of M that is closed under the operation and contains 
the identity element. For instance, the set Nq together with component- 
wise addition and the zero vector forms a commutative monoid and each 
submonoid of it is called a numerical monoid. 

A Hilbert basis of a submonoid K of Nq is a minimal (with respect to 
inclusion) finite subset % oi K such that each element k & K can be written 



(h(/)) = (lt(/) I / G /). 



(2) 



{\t{I)) = {\t{g)\geg). 



(3) 
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as a sum k = Xl/iew '^^h, where Ch € Nq. It is known that each numerical 
submonoid has a unique Hilbert basis [IT] . 

Submonoids arise in various fields like integer programming. Such a prob- 
lem is usually expressed in standard form: 

Minimize c"^x such that Ax. = b, x > 0, (4) 

where b G Z*", c G Z" and A G Z™^" are given and a non-negative integer 
vector X is to be found. The set of all integer vectors x > satisfying the 
constraint equation Ax = b is called the feasible region. Of interest here is 
the case b = because then the feasible region is the kernel of the matrix A, 
written ker(y4), which is clearly a numerical submonoid. The problem is then 
to find a Hilbert basis of the submonoid K = ker(y4) in Nq , where A = (aij) 
is an m X n integer matrix. 

Following [7] we present an algorithm that solves this problem by using 
Grobner bases. This procedure can also be found in [TOl [13] . 

The first step is to translate this problem into the realm of polynomials. 
To this end, we associate a variable Xi to every row of A, 1 < i < m. Since 
entries of A can be negative integers, we have to consider the ring of Laurent 
polynomials. Furthermore, define the mapping 

Ip : K[vi,. . . ,Vn,Wi,. . . ,Wn] ]K[xf \ . . . , x^^] [wi , . . . , w„] (5) 

on the variables first 

m 

ip{vj) = Wj Y\_^T^ = ^i? 1 < j < "^5 (6) 

1=1 

and then extend it such that it becomes a ring homomorphism. In view of 
the ideal 

lA = (w,f[xT^ -v,\l<j<n'j (7) 

in K[xf^, x^^] [vi, . . . ,Vn,wi, . . . , Wn], we have by [3] 

ker{ijj) = IaHKIvi, . . . ,Vn,Wi, . . . ,Wn\. (8) 

Using this notation and the polynomial ring in (JT]) instead of the ring of 
Laurent polynomials, we obtain the following assertion due to [T3] : 
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Let Q he & Grobner basis for with respect to any monomial 
order for which Xi >~ Vj, t >- Vj and Vj >- Wi for all 1 < i < m and 
1 < J < 'T-- A Hilbert basis for K = ker(74) is then given by 

7/ = {a e N[J I v" - w" G 6?} . (9) 

A proof can be found in [13j . 

This result facilitates an algorithm for computing the Hilbert basis of a 
given submonoid ker(y4), which is summarized by Algorithm [TJ 

Algorithm 1 Grobner basis algorithm for computing a Hilbert basis. 

1. Associate the ideal I a defined in ([7]) to a given mxn integer matrix A. 

2. Compute the reduced Grobner basis Q for I a with respect to a mono- 
mial order with Xi >~ vj, t >- vj and vj >- Wk for all 1 < z < m and 
1 < j, k < n. 

3. Read off the elements of the shape v° — w", a G Nq, which form a 
Hilbert basis for ker(y4). 



2.3 Linear Codes 

Let F be the finite field. A linear code C of length n and dimension k over F is 
the image of a one-to-one linear mapping : ^ F", i.e., C = ^(F'^), where 
k<n. The code C is an [n, k] code and its elements are called codewords. In 
algebraic coding, the codewords are always written as row vectors. 

A generator matrix for an [n, k] code C is a /c x n matrix G whose rows 
form a basis of C, i.e., C = {slG \ a G F'^}. The code C is in standard form 
if it has a generator matrix in reduced echelon form G = [I^ \ M), where Ik 
is the k X k identity matrix. Each linear code is equivalent (by a monomial 
transformation) to a linear code in standard form. 

For an [n, k] code C over F, the dual code C"*" is given by all words u G F"' 
such that (u, c) = for each c G C, where (■, ■) denotes the ordinary inner 
product. The dual code C"*- is an [n, n—k] code. If G = (7^ | M) is a generator 
matrix for C, then H = (— | In-k) is a generator matrix for C^. For each 
word c G F", c G C if and only if cH^ = 0. The matrix if is a parity check 
matrix for C P, [T3] . 
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3 A Grobner basis algorithm for finding a 
Hilbert basis of a matrix kernel 



In the following, let ¥p denote a finite field with p elements, where p is prime. 
We are interested in finding the Hilbert basis of the submonoid 

K = keT{Hp)n¥;, (10) 

where H is an m x n integer matrix and Hp = H ®i Fp. 

In other words, we are considering the case in which the numerical monoid 
Nq is replaced by the vector space FJJ over the finite prime field Fp. Then the 
submonoid K becomes a vector space and the Hilbert basis equals an ordinary 
basis in the sense of linear algebra. Clearly, the uniqueness property does no 
longer hold. Nevertheless, the Grobner basis algorithm for finding a Hilbert 
basis as described in the previous section (see Algorithm [T]) can be adapted 
to this situation in order to find one vector space basis. 

Since p is congruent in Fp, the following additional ideal will be used 

Jp(x) = (xf - 1 I 1 < z < n) . 

In this way, the exponents of the monomials can be treated as vectors in Fp. 
Let H = {hij) be an m X n-matrix with entries in Fp and define the ideals 

Jh = (vj - Wj f[ 4- I 1 < J < (11) 

and 



Jj^ + Jp(x) + Jp(v) + Jp(w). (12) 



The homomorphism ip defined in ([5]) and ([6]) can be used to detect ele- 
ments in the kernel of H. However, all entries of H can be written (mod- 
ulo p) as integers in {0, 1, . . . , p — 1} and so the Laurent polynomials become 
ordinary polynomials. Hence, the image of ip lies in the polynomial ring 
K[xi, ... ,Xm][wi, ... ,Wn]- Note that each non-zero vector a G Fp can be 
written as 

a = (0, . . . , 0, a,, a), (13) 
where aj G Fp \ {0} and a G F^"*. Furthermore, put 

a' = ajCj — a = (0, . . . , 0, 0, —a), (14) 
where e,- is the ith unit vector. 
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Lemma 1. Let H be an m x n-matrix with entries in ¥p. For each 
non-zero element a G F^, we have 

a e ker(if) ^ -v"'w") = mod [/p(x) + /p(v) + /p(w)] . 

Proof. All computations will be performed modulo /p(x) +/p(v) + /p(w). 
By the definition of ip, we have 



fc=l 1=1 k=l 



\k=l i=l k=l 

In the second equation, w"'w° = w°'"^" = w*^*"' = w'^'. Thus 

- v°'w°) = ^ x^^'"' - x^°' = 

HeiUi — i/a' = (ejCtj — a') = Ha = 
<^=^ a e ker(iJ). 

□ 

Note that ker(?/') is a toric ideal [3], which can be written as 

ker(V') = JHnK[v,w]. (15) 

Inspired by the assertion on Hilbert bases for numerical submonoids and 
based on the previous lemma, we obtain the following main result. 

Theorem 2. Let Q be a Grobner basis for defined as in / fT2j) with 
respect to the lexicographical order with xi y . . . >- Xm >- Vi >~ . . . >- Vn >- 
Wi y . . . y Wn- Then a basis for ker(if) in is given by 



n 



|(0,...,0,ai,a) G F;J | - v°'w" G ^, a' = a^e^ - a, ^ o| 



Using this assertion, we can obtain an adapted version of Algorithm [T] for 
computing a basis for ker(i7) as a subspace of (see Algorithm [2]) . For the 
proof of correctness, which comes hand in hand with the proof of Theorem 
[21 three facts will be required: 
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Algorithm 2 Grobner basis algoritlim for computing a basis for ker(/J). 



1. Associate the ideal Ih defined as in (1121) to a given m x n-matrix H 
over Fp. 

2. Compute the reduced Grobner basis Q for Ih with respect to the lexi- 
cographical order with xi y . . .yxm >~vi ^ . . . >- 1>„ >~wi y . . .>~Wn- 

3. Read off the elements of the form v'^" — v"'w" with a' = Ojej — a and 
ttj 7^ 0, which give a basis for ker(if). 



1. The reduced Grobner basis of a binomial ideal consists of binomials [8]. 

2. The ideal Jh contains no monomials. 

3. The ideal Jh is prime and Ih resembles a prime ideal in the following 
sense: If f,g E /c[x, v,w] are polynomials such that each variable Xi 
involved in fg has an exponent of at most p — 1, i.e., the exponents 
of the monomials are written as elements in F^, then fg G Ih implies 
either f E Ih or g E Ih- 

The following proof is an adapted version of the one in [13] . Note that all 
subsequently performed calculations will be either in Fp or modulo the ideal 
/p(x) + Jp(v) + /p(w). 

Proof. We need to show that the obtained set is a minimal spanning 
set. Assume that this is not the case. Then there must be a non-zero element 
P G ker(if) that cannot be written as a linear combination of elements in 
"H. Choose an element /3 such that the monomial x'^ is minimal with respect 
to the chosen monomial order. Write /3 = (0. . . . , 0, f3i, (3), where A 7^ and 
/3 G ¥p~\ By Lemma [H (fT5|) . and keT{ip) C Jh, we obtain 

/ = - v^'w^ G Jh- 

Thus / can be reduced to zero on division by Q, since Jh C Ih- Hence 
by the definition of Grobner bases, there must be a polynomial g E Q with 
\i{g) = v]" and 1 < 7^ < f3i. Put S = Pi—'ji. In view of the chosen elimination 
order and the fact that Q consists of binomials, it follows that g is of the form 

g = vf - v^'w''. 
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for some 7' = (0, . . . , 0, —7), where 7 G and f] G F^. But by Lemma [H 

the Grobner basis element g will vanish under ip and so 



7] = 'JiBi + 7' 



7- 



Then we have 



^(0...0 5-7)^{0...0 7i7) 
^(0...0 0-7)^(0 ...0 7,7) 



^(0... 0-/3)^(0 ...Oft /3) 
^5 _ ^(0 ... -,3+7)^(0 ... 5 /3-7) 



^-/3'+7'^/3'-7'+5e, 



Applying the previous stated facts 2 and 3 yields 

_ ^-/3'+7V/^'-V+^e, ^ 

Thus by Lemma [H — 7' + 5ei G ker(if). But by the choice oi g, 5 < f3i and 



so x^'^''''^'^'^' -< x^. Hence by the selection of /9, /3' — 7' + 5ej can be written 
as a linear combination of elements in T-L. The same holds for 7, since it lies 
in T-L due to the choice of the corresponding Grobner basis element g. But 
then 

P = /3' + Ae. = I3'+{S + 7^) e, + 7' - 7' = (/3' - 7' + ^e,) + 7, 

and so (3 can also be written as such a linear combination contradicting the 
choice of (3 and hence proving the assertion. □ 

We conclude by giving an example illustrating applications to linear 
codes. 

Example 1. Consider the [11,6] ternary Golay code [9l [H] with the gener- 
ator matrix G = M) , where 

/I 1 1 1 1\ 
12 2 1 
10 12 2 
2 10 12 
2 2 10 1 
\1 2 2 1 0/ 



M 
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Then a parity check matrix is 

/2 

2 
2 



H 




2 
1 
1 
2 



2 
1 
1 
2 





1 










1 / 



Applying Algorithm [2] for computing a basis of ker(if) yields the following 
polynomials belonging to the reduced Grobner basis 



vrVsVgvf^w^WjwlwgWn, 



V6 

- v'^VgVioVuW-iWYWgwjQwj^, 
V2 - vlvgViovf^W2W8wlwl(^Wn, 
Vl - V^vlvgvlQvliWiWjWsWgWioWn- 

The Hilbert basis taken from these polynomials corresponds to the row vec- 
tors of the matrix G. 
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